Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles.

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles.

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

16S rRNA gene-based meta-analysis of the reptile gut microbiota reveals environmental effects, host influences and a limited core microbiota.

An animal's gut microbiota plays an important role in host health, reproduction and digestion. However, many studies focus on only a few individuals or a single species, limiting our ability to recognize emergent patterns across a wider taxonomic grouping. Here, we compiled and reanalysed published 16S rRNA gene sequence data for 745 gut microbiota samples from 91 reptile species using a uniform bioinformatics pipeline to draw broader conclusions about the taxonomy of the reptile gut microbiota and the forces shaping it. Our meta-analysis revealed the significant differences in alpha- and beta-diversity across host order, environment, diet, habitat and conservation status, with host diet and order contributing the most to these differences. We identified the principal bacterial phyla present in the reptile gut microbiota as Bacteroidota, Proteobacteria (mostly Gamma class), and Firmicutes, and detected the bacterial genus Bacteroides in most reptile individuals, thus representing a putative 'core' microbiota. Our study provides novel insights into key drivers of the reptile gut microbiota, highlights existing knowledge gaps and lays the groundwork for future research on these fascinating hosts and their associated microbes.

Exploring the relationship between environmental DNA concentration and biomass in Asian giant softshell turtle (Pelochelys cantorii).

In recent years, environmental DNA (eDNA) technology has become an accepted approach for investigating rare and endangered species because of its economic efficiency, high sensitivity, and non-invasiveness. The Asian giant softshell turtle (Pelochelys cantorii) is a first-class protected aquatic animal in China, and traditional resource survey methods have not identified its natural populations for many years. In this study, primers and a TaqMan probe targeting ND5 were designed, reaction conditions were optimized, a standard curve was constructed using synthetic DNA, and an eDNA quantitative PCR (qPCR) detection method was established. The eDNA detection technology for P. cantorii revealed that the number of species in the experimental pools showed a significant linear relationship with the eDNA concentration (p < 0.05). The eDNA concentration was negatively correlated with the length of time after the removal of P. cantorii and retention in the water body for 9 days. The qPCR detection method for P. cantorii eDNA established in this study can be applied to the qualitative detection of P. cantorii in water bodies, as well as to preliminary evaluation of its relative biomass. This can serve as a baseline for the investigation of natural P. cantorii population and the evaluation of its wild release effects.

Cysteine Enrichment Mediates Co-Option of Uricase in Reptilian Skin and Transition to Uricotelism.

Uric acid is the main means of nitrogen excretion in uricotelic vertebrates (birds and reptiles) and the end product of purine catabolism in humans and a few other mammals. While uricase is inactivated in mammals unable to degrade urate, the presence of orthologous genes without inactivating mutations in avian and reptilian genomes is unexplained. Here we show that the Gallus gallus gene we name cysteine-rich urate oxidase (CRUOX) encodes a functional protein representing a unique case of cysteine enrichment in the evolution of vertebrate orthologous genes. CRUOX retains the ability to catalyze urate oxidation to hydrogen peroxide and 5-hydroxyisourate (HIU), albeit with a 100-fold reduced efficiency. However, differently from all uricases hitherto characterized, it can also facilitate urate regeneration from HIU, a catalytic property that we propose depends on its enrichment in cysteine residues. X-ray structural analysis highlights differences in the active site compared to known orthologs and suggests a mechanism for cysteine-mediated self-aggregation under H2O2-oxidative conditions. Cysteine enrichment was concurrent with the transition to uricotelism and a shift in gene expression from the liver to the skin where CRUOX is co-expressed with ß-keratins. Therefore, the loss of urate degradation in amniotes has followed opposite evolutionary trajectories: while uricase has been eliminated by pseudogenization in some mammals, it has been repurposed as a redox-sensitive enzyme in the reptilian skin.

An Enhanced Method for the Use of Reptile Skin Sheds as a High-Quality DNA Source for Genome Sequencing.

With the emergence of high-throughput sequencing technology, a number of non-avian reptile species have been sequenced at the genome scale, shedding light on various scientific inquiries related to reptile ecology and evolution. However, the routine requirement of tissue or blood samples for genome sequencing often poses challenges in many elusive reptiles, hence limiting the application of high-throughput sequencing technologies to reptile studies. An alternative reptilian DNA resource suitable for genome sequencing is in urgent need. Here, we used the corn snake (Pantherophis guttatus) as a reptile model species to demonstrate that the shed skin is a high-quality DNA source for genome sequencing. Skin sheds provide a noninvasive type of sample that can be easily collected without restraining or harming the animal. Our findings suggest that shed skin from corn snakes yields DNA of sufficient quantity and quality that are comparable to tissue DNA extracts. Genome sequencing data analysis revealed that shed skin DNA is subject to bacteria contamination at variable levels, which is a major issue related to shed skin DNA and may be addressed by a modified DNA extraction method through introduction of a 30 min pre-digestion step. This study provides an enhanced method for the use of reptile shed skins as a high-quality DNA source for whole genome sequencing. Utilizing shed skin DNA enables researchers to overcome the limitations generally associated with obtaining traditional tissue or blood samples and promises to facilitate the application of genome sequencing in reptilian research.

Gene expression of male pathway genes sox9 and amh during early sex differentiation in a reptile departs from the classical amniote model.

BACKGROUND: Sex determination is the process whereby the bipotential embryonic gonads become committed to differentiate into testes or ovaries. In genetic sex determination (GSD), the sex determining trigger is encoded by a gene on the sex chromosomes, which activates a network of downstream genes; in mammals these include SOX9, AMH and DMRT1 in the male pathway, and FOXL2 in the female pathway. Although mammalian and avian GSD systems have been well studied, few data are available for reptilian GSD systems. RESULTS: We conducted an unbiased transcriptome-wide analysis of gonad development throughout differentiation in central bearded dragon (Pogona vitticeps) embryos with GSD. We found that sex differentiation of transcriptomic profiles occurs at a very early stage, before the gonad consolidates as a body distinct from the gonad-kidney complex. The male pathway genes dmrt1 and amh and the female pathway gene foxl2 play a key role in early sex differentiation in P. vitticeps, but the central player of the mammalian male trajectory, sox9, is not differentially expressed in P. vitticeps at the bipotential stage. The most striking difference from GSD systems of other amniotes is the high expression of the male pathway genes amh and sox9 in female gonads during development. We propose that a default male trajectory progresses if not repressed by a W-linked dominant gene that tips the balance of gene expression towards the female trajectory. Further, weighted gene expression correlation network analysis revealed novel candidates for male and female sex differentiation. CONCLUSION: Our data reveal that interpretation of putative mechanisms of GSD in reptiles cannot solely depend on lessons drawn from mammals.

In vitro culture of reptile PGCS to preserve endangered species.

Primordial germ cells (PGCs), are the source of gametes in vertebrates. There are similarities in the development of PGCs of reptiles with avian and mammalian species PGCs development. PGCs culture has been performed for avian and mammalian species but there is no report for reptilian PGCs culture. In vitro culture of PGCs is needed to produce transgenic animals, preservation of endangered animals and for studies on cell behaviour and research on fertility. Reptiles are traded as exotic pets and a source of food and they are valuable for their skin and they are useful as model for medical research. Transgenic reptile has been suggested to be useful for pet industry and medical research. In this research different aspects of PGCs development was compared in three main classes of vertebrates including mammalian, avian and reptilian species. It is proposed that a discussion on similarities between reptilian PGCs development with avian and mammalian species helps to find clues for studies of reptilian PGCs development details and finding an efficient protocol for in vitro culture of reptilian PG.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates.

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Molecular evolution of the hemoglobin gene family across vertebrates.

Adaptation to various altitudes and oxygen levels is a major aspect of vertebrate evolution. Hemoglobin is an erythrocyte protein belonging to the globin superfamily, and the α-, ß-globin genes of jawed vertebrates encode tetrameric ((α2ß2) hemoglobin, which contributes to aerobic metabolism by delivering oxygen from the respiratory exchange surfaces into cells. However, there are various gaps in knowledge regarding hemoglobin gene evolution, including patterns in cartilaginous fish and the roles of gene conversion in various taxa. Hence, we evaluated the evolutionary history of the vertebrate hemoglobin gene family by analyses of 97 species representing all classes of vertebrates. By genome-wide analyses, we extracted 879 hemoglobin sequences. Members of the hemoglobin gene family were conserved in birds and reptiles but variable in mammals, amphibians, and teleosts. Gene motifs, structures, and synteny were relatively well-conserved among vertebrates. Our results revealed that purifying selection contributed substantially to the evolution of all vertebrate hemoglobin genes, with mean dN/dS (ω) values ranging from 0.057 in teleosts to 0.359 in reptiles. In general, after the fish-specific genome duplication, the teleost hemoglobin genes showed variation in rates of evolution, and the ß-globin genes showed relatively high ω values after a gene transposition event in amniotes. We also observed that the frequency of gene conversion was high in amniotes, with fewer hemoglobin genes and higher rates of evolution. Collectively, our findings provide detail insight into complex evolutionary processes shaping the vertebrate hemoglobin gene family, involving gene duplication, gene loss, purifying selection, and gene conversion.

Evolution of the germline mutation rate across vertebrates.

The germline mutation rate determines the pace of genome evolution and is an evolving parameter itself1. However, little is known about what determines its evolution, as most studies of mutation rates have focused on single species with different methodologies2. Here we quantify germline mutation rates across vertebrates by sequencing and comparing the high-coverage genomes of 151 parent-offspring trios from 68 species of mammals, fishes, birds and reptiles. We show that the per-generation mutation rate varies among species by a factor of 40, with mutation rates being higher for males than for females in mammals and birds, but not in reptiles and fishes. The generation time, age at maturity and species-level fecundity are the key life-history traits affecting this variation among species. Furthermore, species with higher long-term effective population sizes tend to have lower mutation rates per generation, providing support for the drift barrier hypothesis3. The exceptionally high yearly mutation rates of domesticated animals, which have been continually selected on fecundity traits including shorter generation times, further support the importance of generation time in the evolution of mutation rates. Overall, our comparative analysis of pedigree-based mutation rates provides ecological insights on the mutation rate evolution in vertebrates.

The conservation and diversity of the exons encoding the glycine and arginine rich domain of the fibrillarin gene in vertebrates, with special focus on reptiles and birds.

The nucleolar rRNA 2'-O-methyltransferase fibrillarin (FBL) contains a highly conserved methyltransferase domain at the C-terminus and a diverse glycine arginine-rich (GAR) domain at the N-terminus in eukaryotes. We found that a nine-exon configuration of fbl and exon 2-3 encoded GAR domain are conserved and specific in vertebrates. All internal exons except exon 2 and 3 are of the same lengths in different vertebrate lineages. The lengths of exon 2 and 3 vary in different vertebrate species but the ones with longer exon 2 usually have shorter exon 3 complementarily, limiting lengths of the GAR domain within a certain range. In tetrapods except for reptiles, exon 2 appears to be longer than exon 3. We specifically analyzed different lineages of reptiles for their GAR sequences and exon lengths. The lengths of exon 2 in reptiles are around 80-130-nt shorter and the lengths of exon 3 in reptiles are around 50-90 nt longer than those in other tetrapods, all in the GAR-coding regions. An FSPR sequence is present at the beginning of the GAR domain encoded by exon 2 in all vertebrates, and a specific FXSP/G element (X can be K, R, Q, N, and H) exist in the middle of GAR with phenylalanine as the 3rd exon 3-encoded amino acid residue starting from jawfish. Snakes, turtles, and songbirds contain shorter exon 2 compared with lizards, indicating continuous deletions in exon 2 and insertions/duplications in exon 3 in these lineages. Specifically, we confirmed the presence the fbl gene in chicken and validated the RNA expression. Our analyses of the GAR-encoding exons of fbl in vertebrates and reptiles should provide the basis for further evolutionary analyses of more GAR domain encoding proteins.

Na+/Cl- cotransporter 2 is not fish-specific and is widely found in amphibians, non-avian reptiles, and select mammals.

Solute carrier 12 (Slc12) is a family of electroneutral cation-coupled chloride (Cl-) cotransporters. Na+/K+/2Cl- (Nkcc) and Na+/Cl- cotransporters (Ncc) belong to the Nkcc/Ncc subfamily. Human and mouse possess one gene for the Na+/Cl- cotransporter (ncc gene: slc12a3), whereas teleost fishes possess multiple ncc genes, slc12a3 (ncc1) and slc12a10 (ncc2), in addition to their species-specific paralogs. Amphibians and squamates have two ncc genes: slc12a3 (ncc1) and ncc3. However, the evolutionary relationship between slc12a10 and ncc3 remains unresolved, and the presence of slc12a10 (ncc2) in mammals has not been clarified. Synteny and phylogenetic analyses of vertebrate genome databases showed that ncc3 is the ortholog of slc12a10, and slc12a10 is present in most ray-finned fishes, coelacanths, amphibians, reptiles, and a few mammals (e.g., platypus and horse) but pseudogenized or deleted in birds, most mammals, and some ray-finned fishes (pufferfishes). This shows that slc12a10 is widely present among bony vertebrates and pseudogenized or deleted independently in multiple lineages. Notably, as compared with some fish that show varied slc12a10 tissue expression profile, spotted gar, African clawed frog, red-eared slider turtle, and horse express slc12a10 in the ovaries or premature gonads. In horse tissues, an unexpectedly large number of splicing variants for Slc12a10 have been cloned, many of which encode truncated forms of Slc12a10, suggesting that the functional constraints of horse slc12a10 are weakened, which may be in the process of becoming a pseudogene. Our results elaborate on the evolution of Nkcc/Ncc subfamily of Slc12 in vertebrates.NEW & NOTEWORTHY slc12a10 is not a fish-specific gene and is present in a few mammals (e.g., platypus and horse), non-avian reptiles, amphibians, but was pseudogenized or deleted in most mammals (e.g., human, mouse, cat, cow, and rhinoceros), birds, and some ray-finned fishes (pufferfishes).

Comparative analysis of codon usage of mitochondrial genomes provides evolutionary insights into reptiles.

Current available information on reptile genomes provides great potential for the study of unique adaptations from a genomic perspective. We compared differences in base composition and codon usage patterns across 400 reptile mitochondrial genomes assessing AT and GC skew, GC frequency, codon usage, effective number of codons, and codon adaptation index. We identified poor GC content in reptile mitochondrial genomes, with a predominant bias toward Adenine. We determined a compositional asymmetry between different taxonomic groups, which are inversely correlated to the rates of rearrangements in the mitogenome. We found that the most common codons in reptile mitochondrion are CTA (L), ATA (M) and ACA (T), which relates with have been found in birds, meaning that these patterns are shared across sauropsid mitogenomes. Codon usage bias clustering and effective codon number analyses revelated compositional asymmetry based on RSCU as well as that reptile mitogenomes are translationally efficient and are under selection pressure. Codon adaptation index revealed highest values in turtles indicating higher translational efficiency of mitochondrial genes among all reptiles, which could be related to metabolic adaptations (i.e., tolerance to anoxic conditions). This was also seen in other groups such as crocodiles (i.e., acclimation to cold) and snakes (phylogenetic origin of toxin-secreting oral glands and the evolutionary redesign of cytochrome c oxidase complex genes). We discuss our findings in the context of potential adaptations and evolutionary implications that these genomic differences provide to reptiles.

Phylogenetically under-dispersed gut microbiomes are not correlated with host genomic heterozygosity in a genetically diverse reptile community.

While key elements of fitness in vertebrate animals are impacted by their microbiomes, the host genetic characteristics that factor into microbiome composition are not fully understood. Here, we correlate host genomic heterozygosity and gut microbiome phylogenetic diversity across a community of reptiles in southwestern New Mexico to test hypotheses about the behaviour of host genes that drive microbiome assembly. We find that microbiome communities are phylogenetically under-dispersed relative to random expectations, and that host heterozygosity is not correlated with microbiome diversity. Our analyses reinforce results from functional genomic work that identify conserved host immune and nonimmune genes as key players in microbiome assembly, rather than gene families that rely on heterozygosity for their function.

Sex Determination and Ovarian Development in Reptiles and Amphibians: From Genetic Pathways to Environmental Influences.

BACKGROUND: Reptiles and amphibians provide untapped potential for discovering how a diversity of genetic pathways and environmental conditions are incorporated into developmental processes that can lead to similar functional outcomes. These groups display a multitude of reproductive strategies, and whereas many attributes are conserved within groups and even across vertebrates, several aspects of sexual development show considerable variation. SUMMARY: In this review, we focus our attention on the development of the reptilian and amphibian ovary. First, we review and describe the events leading to ovarian development, including sex determination and ovarian maturation, through a comparative lens. We then describe how these events are influenced by environmental factors, focusing on temperature and exposure to anthropogenic chemicals. Lastly, we identify critical knowledge gaps and future research directions that will be crucial to moving forward in our understanding of ovarian development and the influences of the environment in reptiles and amphibians. KEY MESSAGES: Reptiles and amphibians provide excellent models for understanding the diversity of sex determination strategies and reproductive development. However, a greater understanding of the basic biology of these systems is necessary for deciphering the adaptive and potentially disruptive implications of embryo-by-environment interactions in a rapidly changing world.

Heterogeneous Evolution of Sex Chromosomes in the Torrent Frog Genus Amolops.

In sharp contrast to birds and mammals, in numerous cold-blooded vertebrates, sex chromosomes have been described as homomorphic. This sex chromosome homomorphy has been suggested to result from the high turnovers often observed across deeply diverged clades. However, little is known about the tempo and mode of sex chromosome evolution among the most closely related species. Here, we examined the evolution of sex chromosome among nine species of the torrent frog genus Amolops. We analyzed male and female GBS and RAD-seq from 182 individuals and performed PCR verification for 176 individuals. We identified signatures of sex chromosomes involving two pairs of chromosomes. We found that sex-chromosome homomorphy results from both turnover and X-Y recombination in the Amolops species, which simultaneously exhibits heterogeneous evolution on homologous and non-homologous sex chromosomes. A low turnover rate of non-homologous sex chromosomes exists in these torrent frogs. The ongoing X-Y recombination in homologous sex chromosomes will act as an indispensable force in preventing sex chromosomes from differentiating.

Cloacal virome of an ancient host lineage - The tuatara (Sphenodon punctatus) - Reveals abundant and diverse diet-related viruses.

Tuatara (Sphenodon punctatus) are one of the most phylogenetically isolated vertebrate species and provide a unique host system to study virus evolution. While the tuatara genome, sequenced in 2020, revealed many endogenous viral elements, we know little of the exogenous viruses that infect tuatara. We performed a metatranscriptomics study of tuatara cloaca samples from a wild population on Takapourewa (Stephens Island), Aotearoa New Zealand. From these data we identified 49 potentially novel viral species that spanned 19 RNA viral families and/or orders, the vast majority (48) of which were likely dietary-related. Notably, using a protein structure homology search, we identified a highly divergent novel virus within the Picornaviridae which may directly infect tuatara. Additionally, two endogenous tuatara adintoviruses were characterised that exhibited long-term viral-host co-divergence. Overall, our results indicate that the tuatara cloacal virome is highly diverse, likely due to a large number of dietary-related viruses.

Computational analysis on two putative mitochondrial protein-coding genes from the Emydura subglobosa genome: A functional annotation approach.

Rapid advancements in automated genomic technologies have uncovered many unique findings about the turtle genome and its associated features including olfactory gene expansions and duplications of toll-like receptors. However, despite the advent of large-scale sequencing, assembly, and annotation, about 40-50% of genes in eukaryotic genomes are left without functional annotation, severely limiting our knowledge of the biological information of genes. Additionally, these automated processes are prone to errors since draft genomes consist of several disconnected scaffolds whose order is unknown; erroneous draft assemblies may also be contaminated with foreign sequences and propagate to cause errors in annotation. Many of these automated annotations are thus incomplete and inaccurate, highlighting the need for functional annotation to link gene sequences to biological identity. In this study, we have functionally annotated two genes of the red-bellied short-neck turtle (Emydura subglobosa), a member of the relatively understudied pleurodire lineage of turtles. We improved upon initial ab initio gene predictions through homology-based evidence and generated refined consensus gene models. Through functional, localization, and structural analyses of the predicted proteins, we discovered conserved putative genes encoding mitochondrial proteins that play a role in C21-steroid hormone biosynthetic processes and fatty acid catabolism-both of which are distantly related by the tricarboxylic acid (TCA) cycle and share similar metabolic pathways. Overall, these findings further our knowledge about the genetic features underlying turtle physiology, morphology, and longevity, which have important implications for the treatment of human diseases and evolutionary studies.

Energy as the cornerstone of environmentally driven sex allocation.

In recent years, observations of distinct organisms have linked the quality of the environment experienced by a given individual and the sex it will develop. In most described cases, facing relatively harsh conditions resulted in masculinization, while thriving in favorable conditions promoted the development of an ovary. This was shown indistinctively in some species presenting a genetic sex determination (GSD), which were able to sex-reverse, and in species with an environmental sex determination (ESD) system. However, this pattern strongly depends on evolutionary constrains and is detected only when females need more energy for reproduction. Here, I describe the mechanisms involved in this environmentally driven sex allocation (EDSA), which involves two main energy pathways, lipid and carbohydrate metabolism. These pathways act through various enzymes and are not necessarily independent of the previously known transducers of environmental signals in species with ESD: calcium-redox, epigenetic, and stress regulation pathways. Overall, there is evidence of a link between energy level and the sexual fate of individuals of various species, including reptiles, fish, amphibians, insects, and nematodes. As energy pathways are evolutionarily conserved, this knowledge opens new avenues to advance our understanding of the mechanisms that allow animals to adapt their sex according to the local environment.